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1.  Statement  of  the  Problem  Studied 


The  goal  of  the  Army  Research  Office  Geosciences  program  is  to  measure  the  three  dimensional 
wind  field  in  the  planetary  boundary  layer  (PEL)  over  a  measurement  volume  with  a  50  meter 
spatial  resolution  and  with  measurement  accuracies  of  the  order  of  20  cm/sec.  The  objective  of 
this  work  is  to  develop  and  evaluate  a  high  vertical  resolution  lidar  experiment  using  the  edge 
technique  for  high  accuracy  measurement  of  the  atmospheric  wind  field  to  meet  the  ARO 
requirements.  This  experiment  allows  the  powerful  capabilities  of  the  edge  technique  to  be 
quantitatively  evaluated.  In  the  edge  technique,  a  laser  is  located  on  the  steep  slope  of  a  high 
resolution  spectral  filter.  This  produces  large  changes  in  measured  signal  for  small  Doppler  shifts. 
A  differential  fi-equency  technique  renders  the  Doppler  shift  measurement  insensitive  to  both  laser 
and  filter  frequency  jitter  and  drift.  The  measurement  is  also  relatively  insensitive  to  the  laser 
spectral  width  for  widths  less  than  the  width  of  the  edge  filter.  Thus,  the  goal  is  to  develop  a 
system  which  will  yield  a  substantial  improvement  in  the  state  of  the  art  of  wind  profile 
measurement  in  terms  of  both  vertical  resolution  and  accuracy  and  will  provide  a  unique  capability 
for  atmospheric  wind  studies. 


The  wind  field  is  the  single  most  important  atmospheric  state  variable  required  for  understanding 
atmospheric  motion  and  predicting  weather.  High  spatial  and  temporal  resolution  wind 
measurements  from  ground-based  or  airborne  remote  sensing  systems  are  important  in  describing 
the  evolution  of  mesoscale  circulations  associated  with  various  atmospheric  phenomena  including 
deep  and  shallow  convection,  density  currents  and  gravity  wave  activity.  The  data  is  extremely 
important  in  efforts  to  understand  and  model  these  phenomena.  The  accuracy  and  high  vertical 
resolution  of  the  edge  technique  will  also  provide  a  unique  capability  to  directly  observe  vertical 
winds  at  the  10  to  20  cm/sec  level  which  will  allow  studies  of  turbulent  motion  and  the  evolution 
of  convective  cells  in  the  boundary  layer.  Such  measurements  could  also  be  applied  for  detection 
of  wind  shear^  microbursts,  and  turbulence  --  all  of  which  represent  significant  aircraft  safety 
hazards. 
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n.  Summary  of  Results 


A.  Lidar  Development 

A  ground-based  Doppler  lidar  system  was  developed  under  this  program  to  demonstrate  the 
capability  of  the  edge  technique  to  measure  atmospheric  winds.  In  the  first  part  of  the  program, 
the  hardware  and  software  for  acquisition  of  the  wind  lidar  data  was  developed.  This  work 
utilizes  our  prior  work  developing  the  theory  and  methodology  of  the  edge  technique.  A  brief 
overview  of  the  methodology  of  the  measurement  and  a  description  of  the  experiment  is  provided 
below. 

An  optical  schematic  of  the  lidar  system  is  shown  in  Figure  1.  An  injection  seeded  pulsed 
Nd:YAG  laser  is  used  as  the  narrowband  energy  source.  The  laser  has  an  output  energy  of  up  to 
800  mJ  per  pulse  at  a  repetition  rate  of  10  Hz.  The  laser  pulse  length  is  15  nsec  and  we  have 
measured  the  spectral  width  of  this  laser  as  35-40  MHz  at  1.06  pm  which  is  near  transform 
limited  for  a  Gaussian  pulse.  The  spectral  bandwidth  of  this  laser  is  sufficiently  narrow  to  make 
the  proposed  wind  measurements  with  the  edge  technique. 

A  small  fraction  of  the  outgoing  laser  signal  is  picked  off  and  used  to  make  a  reference 
measurement  of  the  laser  frequency  on  the  Fabry  Perot  etalon  for  each  outgoing  laser  pulse.  The 
beam  is  first  collimated  and  the  divergence  is  matched  to  that  of  the  atmospheric  backscattered 
signal  collected  by  the  telescope.  It  is  also  attenuated  to  provide  signals  at  the  detectors  of  the 
same  magnitude  as  those  backscattered  from  the  atmosphere.  The  reference  beam  is  then  input  to 
the  edge  detection  setup. 

The  major  portion  of  the  laser  energy  (99%)  is  transmitted  to  the  atmosphere  via  beam  steering 
optics.  A  portion  of  the  laser  energy  is  backscattered  from  each  range  element  of  the  atmosphere 
and  the  resultant  signal  is  collected  by  a  reflective  telescope  with  a  diameter  of  0.4  meters.  The 
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field  of  view  of  the  telescope  is  set  by  a  field  stop  which  is  adjusted  to  make  the  field  of  view  of 
the  telescope  greater  than  or  equal  to  the  divergence  of  the  outgoing  laser  beam.  The  resultant 
beam  is  collimated  and  transferred  to  a  beamsplitter  which  combines  the  atmospheric 
backscattered  and  the  picked  off  reference  laser  signals.  The  signals  are  then  bandpassed  by  a 
narrowband  interference  filter  to  reject  unwanted  background  and  are  then  collimated  to  satisfy 
the  collimation  and  throughput  requirements  of  the  Fabry  Perot  etalon  which  is  used  as  the  edge 
filter.  The  collimated  beam  is  then  split  between  the  edge  filter  and  energy  monitor  paths. 
Condensing  optics  focus  the  beams  to  minimize  the  size  of  the  high  quantum  efficiency  solid  state 
avalanche  photodiode  detectors. 

The  central  fnnge  of  a  Fabry-Perot  etalon  is  used  as  a  high  resolution  edge  filter  to  measure  the 
shift  of  the  aerosol  return  which  has  the  original  sharp  spectral  distribution  of  the  laser.  Various 
other  filters  with  sharp  spectral  edges  could  be  used.  These  include  gratings,  prisms,  and 
molecular  or  atomic  gas  absorption  lines.  Fabry  Perot  etalons,  plane  or  spherical,  offer  much 
higher  resolution  (sharper  edges)  than  the  other  filters  and  they  also  have  high  throughput.  In 
addition  the  etalon  bandpass  can  be  tuned  to  match  the  laser  frequency  which  allows  existing 
lasers  such  as  Nd:YAG  to  be  used.  Also,  a  filter  with  a  sharp  transmission  feature  such  as  an 
etalon  has  the  advantage  of  rejecting  broadband  background  sources  such  as  the  Rayleigh 
backscatter  or  solar  background  while  an  atomic  or  molecular  absorption  feature  would  transmit 
these  unwanted  background  sources. 

The  edge  and  energy  monitor  detector  signals  are  simultaneously  sampled  as  a  function  of  time  to 
provide  range  resolved  information.  Two  synchronized  12-bit  analog  to  digital  converters  with 
sampling  rates  as  high  as  20  MHz  are  used  to  obtain  range  resolution  as  high  as  7.5  meters 
(round-trip).  The  collected  signals  are  digitized,  stored  and  processed  for  the  Doppler/velocity 
information. 

The  laser  signal  is  transmitted  coaxially  with  the  telescope  by  a  small  flat  to  beam  steering  optics 
located  outside  our  laboratory.  The  beam  steering  optics  provide  the  capability  to  direct  the  beam 
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to  the  atmosphere  and  to  vary  the  azimuth  and  elevation  angles  of  the  beam.  The  azimuth  and 
elevation  angles  are  manually  adjusted.  The  beam  can  be  directed  from  30  to  120  degrees  in 
elevation  angle  and  from  -60  to  60  degrees  in  azimuth.  This  allows  us  to  make  atmospheric  slant 
path  measurements  to  measure  the  horizontal  components  of  the  wind  field  as  well  as  to 
independently  measure  the  vertical  component  of  the  wind. 

B.  Wind  measurements 

Over  the  past  two  years  we  have  built  a  ground-based  lidar  system  as  described.  Recently  we 
were  successful  in  integrating  and  testing  the  Fabry  Perot  etalon  edge  filter  in  the  lidar  system. 
Tests  were  performed  using  both  a  narrowband  pulsed  and  a  CW  Nd:  YAG  laser  and  the  results 
showed  the  desired  high  resolution  etalon  fnnges.  The  etalon  had  a  plate  separation  of  5  cm 
yielding  a  free  spectral  range  of  0. 1  cm'*  and  a  plate  reflectivity  of  93.5  %  giving  a  theoretical 
finesse  of  47  and  a  theoretical  resolution  (FWHH)  of  65  MHz.  The  observed  fnnges  for  the  CW 
laser  (linewidth  <  100  kHz)  had  a  measured  finesse  as  high  as  44.  The  firinges  observed  for  the 
atmospheric  backscattered  aerosol  signal  had  a  finesse  of  30  which  is  in  good  agreement  with 
theory  when  the  pulsed  laser  linewidth  of  40  MHz  and  the  effects  of  finite  angular  field  of  view 
are  included.  A  finesse  of  30  corresponds  to  a  total  instrument  resolution  (FWHH)  of  100  MHz 
giving  a  sensitivity  of  4%  per  m/s  at  the  half-width  of  the  fringe. 

In  order  to  make  wind  measurements  on  the  edge,  a  servo  control  system  was  developed  to  lock 
the  edge  of  the  etalon  fnnge  to  the  pulsed  Nd;YAG  laser  frequency.  The  reference  measurement 
of  the  outgoing  laser  frequency  on  the  edge  is  compared  to  a  setpoint  value,  for  instance  the  half¬ 
width.  A  deviation  from  the  setpoint  value  generates  a  small  error  signal  which  is  applied  to  tune 
the  etalon  to  compensate  for  the  deviation.  The  purpose  of  this  system  is  to  keep  the  laser 
frequency  at  approximately  the  correct  location  on  the  edge.  Very  high  frequency  stability  is  not 
required  since  a  differential  frequency  measurement  is  made.  That  is,  the  frequency  of  the 
atmospheric  backscattered  signal  from  each  range  gate  is  measured  relative  to  the  frequency  of 
the  outgoing  laser  signal,  the  reference  measurement,  for  each  pulse.  The  time  history  of 
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measurements  on  the  edge  at  a  fixed  altitude  of  800  m  is  shovm  in  Figure  2  for  a  10  shot  average 
and  a  range  resolution  of  22.5  m.  The  normalized  signals  (Iedg/Iem)  for  the  atmosphere  and  the 
corresponding  signals  for  the  outgoing  reference  frequency  are  shown  along  with  the  difference 
signal.  As  may  be  seen,  the  reference  measurement  on  the  edge  varies  as  a  function  of  time. 
However,  the  atmospheric  measurements  track  the  reference  measurements  and  by  taking  the 
difference  between  these  the  pulse-to-pulse  frequency  variations  due  to  the  laser  or  the  filter  are 
removed.  The  difference  measurement  is  proportional  to  the  Doppler  shift  and  contains  the  wind 
information. 

Figure  3  shows  our  first  wind  profile  measurements  which  were  made  on  April  26,  1994.  The 
data  are  for  a  100  shot  average  and  a  range  resolution  of  22.5  m  which  yields  a  vertical  resolution 
of  approximately  15  meters.  The  measurements  are  for  the  radial  component  of  the  wind  along 
the  line  of  sight  of  the  laser  in  the  east-west  direction.  The  Nd:YAG  laser  has  an  output  energy  of 
up  to  800  mJ  per  pulse.  The  laser  energy  actually  used  for  the  wind  measurements  was  limited  to 
120  mJ  to  avoid  damaging  the  metallic  coatings  on  the  40  cm  diameter  atmospheric  beam  steering 
mirror.  The  measurements  were  made  by  manually  steering  the  atmospheric  output  mirror  to 
direct  the  beam  at  an  elevation  angle  of  45  degrees  first  toward  the  east  and  then  towards  the 
west.  The  east- west  component  of  the  wind  was  found  by  taking  half  of  the  difference  in  the 
measured  line  of  sight  winds  for  the  two  directions.  The  rms  error  of  the  data  is  approximately 
0.25  m/s  which  corresponds  to  a  peak-to-peak  variation  of  1.25  m/s. 

In  the  last  six  months  of  the  second  year  of  this  program,  we  made  measurements  with  the  Edge 
Technique  Lidar  system  which  demonstrate  some  of  the  unique  capabilities  of  our  system.  Three 
different  types  of  measurements  were  obtained: 

1)  hard  target  measurements  which  provide  zero  velocity  validation  of  the  lidar; 

2)  profiles  of  vector  winds  i.e.  horizontal  wind  speed  and  direction  and  comparison  of  the 
lidar  wind  profiles  with  rawinsonde  measurements  from  Dulles  airport. 

3)  line-of-sight  wind  profiles  as  a  function  of  time  which  provide  a  high  resolution 
view  of  the  temporal  evolution  of  winds  and  turbulence  in  the  boundary  layer. 
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The  hard  target  measurements  were  undertaken  in  order  to  provide  a  zero  velocity  calibration 
measurement  which  is  a  pivotal  point  for  the  lidar  instrument.  This  meaurement  allows  us  to 
determine  if  there  are  any  systematic  instrument  effects  and  provides  a  statistical  measurement  of 
the  lidar  accuracy  for  zero  velocity. 

The  normalized  edge  signals  measured  by  the  lidar  system  are  shown  in  Figure  4  as  a  function  of 
time.  Each  measurement  represents  a  100  shot  average  (10  seconds).  The  data  represented  by 
circles  are  the  reference  measurements  of  the  outgoing  laser  frequency.  These  are  made  by 
locating  the  position  of  the  outgoing  laser  frequency  on  the  edge  of  the  spectral  response  function 
of  the  Fabry-Perot  etalon  edge  filter.  The  triangles  correspond  to  zero  velocity  hard  target 
measurements  obtained  from  trees  located  about  50  meters  from  the  lidar.  The  position  of  the 
laser  on  the  edge  is  held  constant  from  shot  to  shot  by  locking  the  etalon  to  the  laser  frequency 
with  a  servo  control  loop.  This  eliminates  sensitivity  to  frequency  instabilities  of  the  laser  or  the 
etalon.  We  have  demonstrated  locking  the  etalon  to  the  laser  in  the  presence  of  large  frequency 
drifts  (equivalent  to  >  100  m/s  velocity)  for  many  hours  of  continuous  operation.  This  is 
described  in  our  recent  Applied  Optics  paper.  As  shown  in  Figure  4,  the  lidar  hard  target  return 
overlays  the  reference  data  very  closely  and  tracks  its  value  as  the  reference  varies  slightly  from 
measuremnt  to  measurement.  This  is  what  one  would  expect  to  see  for  a  zero  velocity  target  (i.e. 
no  Doppler  shift). 

Figure  5  shows  the  velocity  measurements  resulting  from  the  hard  target  measurements.  The 
Doppler  shift,  and  thus  the  wind,  is  determined  by  taking  a  difference  between  the  target  and 
reference  measurements  (see  Figure  4)  and  dividing  by  the  sensitivity  (slope)  of  the  edge  at  the 
position  of  the  measurement.  The  sensitivity  for  this  experiment  is  3%  per  m/sec.  That  is,  a  1 
m/sec  wind  produces  a  change  in  measured  signal  of  3%  .  The  absolute  accuracy  of  these 
measurements  is  very  high  with  a  mean  offset  of  only  19  cm/s  and  a  rms  standard  deviation  of 
only  17  cm/s.  The  excellent  agreement  of  these  data  with  the  zero  velocity  hard  target  indicate 
that  the  lidar  system  is  performing  as  expected  and  that  systematic  effects  are  small. 
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Figure  6  is  an  example  of  a  wind  profile  obtained  with  the  lidar  system  showing  wind  speed  and 
direction  as  a  function  of  altitude.  The  horizontal  wind  speed  and  direction  are  obtained  by 
analyzing  two  line-of-sight  wind  profile  measurements  taken  at  45®  elevation  angle  with  a  90° 
difference  in  azimuth  angle  i.e.  one  profile  is  taken  pointing  to  the  East  with  a  second  profile 
pointing  to  the  North.  The  data  are  for  a  100  shot  average  (10  seconds)  and  a  vertical  resolution 
of  26  meters  (37.5  m  range  resolution).  We  estimate  from  the  signal-to-noise  in  the  data  that  the 
accuracy  of  these  wind  measurements  is  about  20  cm/s.  The  Dulles  airport  rawinsonde  wind 
profiles  taken  at  1200  Z  and  0000  Z  on  the  same  day  are  also  shown  in  the  figure  for  comparison. 
The  rawinsonde  profiles  show  only  modest  variations  during  the  12  hour  interval  between 
measurements  which  indicates  fdrly  stable  synoptic  scale  conditions  prevailed.  Under  these 
conditions  the  Dulles  measurements  probably  are  representative  of  the  winds  at  Goddard.  Very 
good  agreement  is  observed  between  the  lidar  and  rawinsonde  measured  winds. 

Figure  7  shows  a  three  dimensional  view  of  the  wind  as  a  function  of  range  (altitude)  and  time 
obtained  by  pointing  the  lidar  in  one  direction  (North)  at  a  45  °  elevation  angle  and  accumulating 
profiles  every  10  seconds.  These  data  represent  a  100  shot  (10  second)  average  and  a  vertical 
resolution  of  26  meters.  The  surface  thus  gives  a  picture  of  the  time  history  of  the  line-of-sight 
wind  and  turbulence  in  the  atmosphere  with  extremely  high  spatial  and  temporal  resolution.  Each 
individual  wind  profile  shows  a  variation  of  the  wind  wdth  range  (altitude)  with  several  prominent 
wind  shear  features  clearly  evolving  during  the  observation  time  of  about  5  minutes.  The  rapidly 
evolving  dynamic  nature  of  the  planetary  boundary  layer  is  well  described.  The  figure  clearly 
shows  the  evolution  of  the  wind  field  in  the  temporal  and  spatial  domains. 

The  last  group  of  figures,  Figures  8-16,  show's  the  line-of  sight  wind  data  in  more  detail  by  displaying 
the  individual  profiles  in  groups  of  3  consecutive  measurements,  one  taken  every  10  seconds.  The 
starting  time  of  the  measurements  is  shown  in  the  upper  right  hand  comer.  The  circles  are  wind 
profile  data  which  result  from  the  first  10  second  average,  the  stars  are  the  next  10  second  average 
(20  seconds  elapsed  time)  and  the  triangles  are  the  third  10  second  average  (30  seconds  elapsed 
time).  The  vertical  resolution  is  still  26  meters.  This  more  conventional  2-dimensional  display  of  the 
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data  gives  a  clear  picture  of  the  rapidly  varying  structure  of  the  wind  field.  The  evolution  of  distinct 
shear  layers  is  readily  apparent.  These  features  are  observed  from  one  independent  profile  to  the 
next.  Of  particular  interest  is  the  evolution  of  the  wind  shear  feature  near  850  m  for  times  from  4;49 
to  4:52.  The  various  shear  features  show  slow  and  consistent  variation  which  is  both  scientifically 
interesting  and  indicates  the  validity  of  the  data  since  this  type  of  slowly  varying  structure  along 
independent  profiles  is  not  consistent  with  any  systematic  instrumental  effects.  Figure  7  also  clearly 
shows  the  evolution  of  the  wind  shear  feature  near  850  m. 
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ND:YAG  LASER 


Figure  i  -  Optical  schematic  of  the  edge  technique  hclar  system. 


Figure  2  -  Normalized  edge  signals  for  the  outgoing  laser  frequency  reference  measurement  and  for 
the  atmospheric  return  from  an  altitude  of  800  m.  Each  data  point  is  determined  from  a  10  shot 
average.  The  difference  (Dev)  contains  the  Doppler  shift  information  from  wliich  the  wind  is 


gurc  3  -  East- West  radial  wind  profile  measurements  obtained  with  the  edge  technique  hdar.  The  data  are  for 
vertical  resolution  of  15  m,  a  100  shot  average,  a  45  deg  elevation  angle  and  a  laser  energy  of  120  mj. 
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Figure  6  -  Lidar  profiles  of  wind  speed  and  direction  obtained  with  the  edge  technique  lidar.  Rawinsonde 
profiles  taken  at  Dulles  airport  on  the  same  day  are  shown  for  comparison. 
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Figure  8  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:49:00  p.m. 
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Figure  9  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:49:30  p.m. 
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Figure  10  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:50:00  p.m. 
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Figure  12  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:51:00  p.m. 
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Figure  13  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:51 :30  p.m. 
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Figure  14  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:52:00  p.m. 
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Figure  15  -  Line  of  sight  wind  profiles  obtained  with  the  edge  technique 
lidar.  Three  consecutive  profiles,  each  taken  from  a  10-second  average, 
are  shown  beginning  at  4:52:30  p.m. 


